figure
subplot(4,3,2), plot(aux_data(8,:))
title('Mouse I Day 2 Rewards')
axis normal
box off
ylim([-0.1 5])
subplot(4,3,5), plot(aa)
title('Mouse I Day 2 Running')
box off
subplot(4,3,1), plot(aux_data4(5,:))
title('Mouse IDay 1 Rewards')
box off
ylim([-0.1 5])
subplot(4,3,4), plot(dd)
title('Mouse I Day 1 Running')
box off
subplot(4,3,3), plot(aux_data5(8,:))
title('Mouse I Day 3 Rewards')
box off
ylim([-0.1 5])
subplot(4,3,6), plot(ee)
title('Mouse I Day 3 Running')
box off


subplot(4,3,8), plot(aux_data2(8,:))
title('Mouse II Day 2 Rewards')
box off
ylim([-0.1 5])
subplot(4,3,11), plot(bb)
title('Mouse II Day 2 Running')
box off
subplot(4,3,7), plot(aux_data3(6,:))
title('Mouse II Day 1 Rewards')
box off
ylim([-0.1 5])
subplot(4,3,10), plot(cc)
title('Mouse II Day 1 Running')
box off
subplot(4,3,9), plot(aux_data6(8,:))
title('Mouse II Day 3 Rewards')
box off
ylim([-0.1 5])
subplot(4,3,12), plot(ff)
title('Mouse II Day 3 Running')
box off

aux_data=load_lvd('D:\MHTempBehavior\S1-T9560.lvd');
aux_data2=load_lvd('D:\MHTempBehavior\S1-T9559.lvd');
aux_data3=load_lvd('D:\MHTempBehavior\S1-T9558.lvd');
aux_data4=load_lvd('D:\MHTempBehavior\S1-T9557.lvd');
aux_data5=load_lvd('D:\MHTempBehavior\S1-T9577.lvd');
aux_data6=load_lvd('D:\MHTempBehavior\S1-T9590.lvd');

[velP,velM,velP_ind,velM_ind,velP_raw,velM_raw]=get_vel_ind_from_adata_M1(aux_data);
aa=smooth2(velP,1000);
aa(aa<0)=0;

[velP,velM,velP_ind,velM_ind,velP_raw,velM_raw]=get_vel_ind_from_adata_M1(aux_data2);
bb=smooth2(velP,1000);
bb(bb<0)=0;

[velP,velM,velP_ind,velM_ind,velP_raw,velM_raw]=get_vel_ind_from_adata_M1(aux_data3);
cc=smooth2(velP,1000);
cc(cc<0)=0;

[velP,velM,velP_ind,velM_ind,velP_raw,velM_raw]=get_vel_ind_from_adata_M1(aux_data4);
dd=smooth2(velP,1000);
dd(dd<0)=0;

[velP,velM,velP_ind,velM_ind,velP_raw,velM_raw]=get_vel_ind_from_adata_M1(aux_data5);
ee=smooth2(velP,1000);
ee(ee<0)=0;

[velP,velM,velP_ind,velM_ind,velP_raw,velM_raw]=get_vel_ind_from_adata_M1(aux_data6);
ff=smooth2(velP,1000);
ff(ff<0)=0;
